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Abstract. We present homogeneous age determinations for a large sample of 55 Galactic globular clusters, which 
constitute about 30% of the total Galactic population. A study of their age distribution reveals that all clusters 
from the most metal poor ones up to intermediate metallicities are coeval, whereas at higher [Fe/H] an age spread 
exists, together with an age-metallicity relationship. At the same time, all clusters within a certain galactocentric 
distance appear coeval, whereas an age spread is present further away from the Galactic centre, without any 
correlation with distance. The precise value of [Fe/H] and galactocentric distance for the onset of the age spread 
and the slope of the age-metallicity relationship are strongly affected by the as yet uncertain [Fe/H] scale. We 
discuss how differences in the adopted [Fe/H] scale and cluster sample size may explain discrepant results about the 
clusters age distribution reached by different authors. Taking advantage of the large number of objects included in 
our sample, we also tested the possibility that age is the global second parameter which determines the Horizontal 
Branch morphology, and found indications that age could explain the global behaviour of the second parameter 
effect. 
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1. Introduction 

Less than a decade ago the age of the oldest globular clus- 
ters appeared to be much higher than that of the expand- 
ing universe. But at the end of the last millenium signifi- 
cant improvements both in models and in observational 
data, notably in the determination of cluster distances 
by virtue of Hipparcos-based distances, lead to a reduc- 
tion of cluster ages. Presently, most determinations scat- 
ter around a typical age of the oldest objects of 12-14 Gyr 
(^alar is fc Weiss, 1998|; [Carretta et al., 2000^ [VandenBcrg 



age determinations in four metallicity ranges (see also 
Richer et al. 1996 for a similar technique). It does not 



200C ; Chaboyer, 2001). With the growing confidence in the 
absolute age determinations and an increasing number of 
extensive homogeneous and high-quality photometric clus- 
ter data, the interest has shifted to questions concerning 
relative ages in order to learn about the formation of the 
galaxy and its halo and disk components. 

In our own work ( palaris et al., 1997 ; Salaris & Weiss 



1997, 1998) we have used a method (to be fully described 



in Sect. 2), which is a mixture of absolute and relative 
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need the knowledge of cluster distances, but instead pre- 
dicts them. To determine the absolute ages of a sample 
of reference clusters for the various metallicity ranges, we 
use the brightness difference between horizontal branch 
(HB) and turn off (TO); the HB sets the cluster distances 
and the TO (whose brightness is the best predicted quan- 
tity to be used as age indicator) their age. For relative 
ages with respect to reference clusters we consider the dif- 
ferences in the extension of the subgiant branch, i.e., the 
colour range between the TO and the base of the red gi- 
ant branch (RGB). In this way we are able to determine 
ages also in case of clusters for which the HB brightness 
determination is problematic (e.g., blue and/or sparsely 
populated HB). Since the differential properties of the TO- 
RGB colour difference across the entire cluster metallic- 
ity range is in principle subject to current uncertainties 
in convection treatment and colour transformations, the 
method works best (i.e., the uncertainties are minimized) 
if restricted to small [Fc/H] ranges. When selecting the 
limits of the individual [Fe/H] intervals, one must strike a 
balance between having a sizable number of clusters, and 
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at the same time covering a [Fe/H] interval where the pre- 
dicted variation of the derivative of the TO-RGB colour 
difference with respect to age and [Fe/H] is small. 

In Salaris & Weiss (1998) we applied this method 
to a sample of 31 clusters taken from various sources, 
such that we could start to address the question of age- 
metallicity and age-galactocentric distance relations. This 
cluster sample appeared still too small to draw statisti- 



method. Evidently, semi-empirical methods for relative 
cluster ages spanning the whole metallicity range reveal 
different and internally inconsistent results for different 
isochrone sets, at least when one does not consider the 
HB brightness from ZAHB models consistent with the em- 
ployed isochrones. 

In this paper we will apply all these methods to a 



cally rieaningful conclusions; strictly speaking, relative 



ages are valid only within each metallicity bin, but we 
have shown in Salaris fc Weiss (1998[) and Pulonc ct al 



cluster sample which is a combination of the Rosenberg 
et al. (20004 |b|) set (35 clusters) and 20 additional clusters 



(1998|y that they arc consistent all over the [Fe/H] range 
spanned by the clusters. In the present paper, we want to 
derive ages for a much larger sample of globular clusters 
by using our own isochrones and applying our preferred 
method. Relative ages for these clusters are a straightfor- 
ward by-product. 

There are other age determination techniques in use, 
some of them being easier to apply, but probably less ac- 
curate. All, of course, depend at some s tage on theoreti- 
cal isochrones. Semi-empirical methods ( Buonanno et al. 
199S |; Ifiosenberg et al., 199S| ) are of particular interest 
for relative age-dating; they make extensive use of the 
colour difference between TO and RGB across all the 
[Fe/H] range spanned by the galactic globular clusters, 
thus avoiding the necessity of RR Lyrae stars or a well- 
developed HB for distance determination. This allows to 
apply the method to a much larger number of clusters, and 



with accurate photometry obtained recently. For this, we 
need colour transformations to {V ~ I); in Weiss & Salaris 
(1999D we investigated the properties of our own choice of 



transformations, showing that for selected clusters with 
photometry in V, B, and / we can obtain identical ages 
and similar isochrone fitting quality for VI and BV photo- 
metric data. In the present paper we will determine cluster 
ages from T//-data using this and a second set of transfor- 
mations ( Alonso et al., 1999 ). We will demonstrate that 
the semi-empirical relative ages across the whole cluster 
metallicity range we obtain from the TO-RGB colour dif- 
ference depend strongly on the chosen transformations; 
those by Alonso ct al. (1999) provide results consistent 
with the ages inferred from the TO-HB brightness differ- 
ences and the method employed in 3alaris & Weiss (1998), 
while the transformations used in Weiss & Salaris (1999) 
are not suited for this purpose. 

On the other hand, as will become evident, our pre- 



ferred age determination method, as used by Salaris & 



thus to obtain a solid picture of cluster formation through- Weiss (1998), turns out to be largely insensitive to the 



out the Galactic history. The calibration of the TO-RGB 
colour difference with respect to age and metallicity is 
obtained empirically making use of clusters of different 
ages and [Fe/H], for which reliable ages from the TO-HB 
brightness difference can be determined. Buonanno et al. 
(1998|)]noticed inconsistencies in relative cluster ages de- 
rived from brightness differences (between HB and TO) or 
colour diff erences as predicted by the isochrones they em- 
ployed. In Pulonc et al. (1998 ) we could demonstrate the 
self-consistency of our own isochrones for (B — y)-colours. 

Following on this, Rosenberg et al. (1999 ) applied the 
semi-empirical method for colour differe nces in {V — I) t o 
their own large sample of cluster data ( Rosenberg et al 
2000a|jp|). Their main result, based on the isochrones by 



Straniero et al. (1997 ) and a pre-release of VandcnBcrg 
et al. ( |2000D models, was that clusters of [Fe/H] < -1.2 



specific colour transformation adopted, although it does 
make use of the colour difference between TO and RGB 
within restricted metallicity ranges. Similarly, the use of 
either VI or BV-data does not influence the picture of the 
absolute and relative ages, but may affect only some indi- 
vidual clusters. Quantitatively, the effect is comparable to 
that of using different photometric sources; it is important 
to realize that in spite of the large observational progress 
and effort, the source of photometric data still influences 
individual and global results on cluster ages. 

The paper is organized as follows: In the next section, 
our preferred age determination method and its appli- 
cation to a large sample of 55 globular clusters will be 
presented. This is, to our knowledge, the largest sample 
of clusters ever investigated, using the same age-dating 
method and theoretical isochrones. The sample size is a 



are pr obably coeval (within 1 Gyr) and that at higher 



metalli cities, up to [Fe/H] < —1.0 a few clusters younger 
by ~ 4 Gyr exist, too, and finally, that at [Fe/H] « —0.7 
younger clusters prevail. However, if the brightness differ- 



significant increase over the 43 clusters by Chaboyer et al 
(1996| ); analysis of large clusters sample h ave been re- 



ccntly performed also b y Richer et al. (1996, 36 objects) 



ence b jtween HB (whose absolute magnitude was taken 



from the empirical relation by Carretta et al. 2000) and 
TO was used as the age indicator, this picture was not 
reproduced completely (some small differences in the age 
distribution for the most metal poor and most metal rich 
clusters), and, most importantly, the absolute age cal- 
ibration for both relative age scales differed by up to 



Buonanno et al. (199S, 33 objects) and Rosenberg et al 
(1999| , 35 objects). Our sample is largely based on ho- 



3 Gyr. A further set of isochrones tested ( Cassisi et al 
1998|)~^howed a large age spread for the colour-difference 



mogeneous and very recent photometries; the isochrones 
and colour transformations are of the latest generation 
as well. In Sect. 3 we will then turn to the semi-empirical 
relative-age method using (F — /)-colours. We will present 
completely consistent results that agree well with those of 
Sect. 2. The discussion in Sect. 4 will analyze the age dis- 
tribution of our cluster sample in light of the proposed sce- 
narios of Galaxy formation. We will present also an anal- 
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ysis of the effect of tlie sample size and the adopted clus- 
ters metallicities on the derived age distribution. About 
the latter point, we have considered cluster [Fe/H] val- 
ues derived by both Carretta fc Gratton (1997| our ref- 
erence metallicities) and Zinn fc West (198^ The two 
sets of values show differences up to ~0.3 dex (Carretta 



& Gratton 1997 [Fe/H] being higher) for intermediate 
metallicity clusters; a qua dratic relationship ([Carretta fc 



Gratto |n 1997[ ) transforms |Zinn fc West (1984|) metaUic- 
ities into the corresponding Carretta fc Gratton (1997) 
ones. The possibility that age is the global second param- 
eter responsible (together with [Fe/H]) of the cluster HB 
morphology will be also tested. A short summary of the 
main results follows in Sect. 5. 

2. The cluster sample and age-dating method 

Our GC sample consists of 55 objects. The V magnitude 
difference between Zero Age HB (ZAHB) and TO, AV, 
and the {V — I) colour difference {VI Johnson-Cousins 
bands) betw een TO and the base o f the RGB, A{V - I) 
as defined in pR.osenberg ct al. (1999 , R99), for 35 of them 
are taken from the homogeneous database by R99. For 
the remaining 20 objects AV^, A{V - I) or A{B - V), 
the corresponding colour difference between TO and base 



of the RGB in {B — V) as defined by VandenBerg et al 



|(1990| ) j are taken from [Safaris fc Weiss (1998|, SW98) and 
other papers published in the last 5 years (see Table |1[) . 

The age-dating method is analogous to the procedure 
outfined in ^alaris fc Weiss (1997| , SW97) and SW98. 
The GC sample is divided into four metallicity inter- 
vals, which are the same ones as in SW98, apart from 
the last interval; this we extended to include in the 
same group the thick disk clusters present in the sample. 
The metallicity intervals therefore are —2.15 <[Fe/H]< 
-1.75, -1.74 <[Fe/H]< -1.3, -1.29 <[Fe/H]< -0.9, and 
-0.89 <[Fe/H]< -0.6 . We use d the [Fe/H] values given by 
Rutledge et al. (1997[) on the [Carretta fc Gratton (1997] , 
CG97) scale (the internal accuracy of these [Fe/H] val- 
ues is of the order of 0.10 dex); this introduces small dif- 
ferences with respect to the [Fe/H] values used by SW98 
whi ch were either direct ly from CG97 or transformed from 
the ^inn fc West (1984[, ZW84) scale using the relation- 



M 3 have replaced M 68 (NGC 4590) and NGC 6584 as 
reference clusters in their respective groups (SW98). Due 
to the small [Fe/H] differences with respect to SW98, the 
cluster NGC 6171 has now moved to the third metallicity 
group where it has replaced M 5 (NGC 5904) as the ref- 
erence cluster. For all the four reference clusters we used 
AV values taken from R99. 

Once the absolute ages of the reference clusters were 
obtained, the age of all other clusters within each group 
was determined differentially with respect to the reference 
one, by using the differences in A{V — I) to derive differ- 
ences in age, again by using the isochrones. In case of clus- 
ters with BV photometry, their relative age with respect 
to the template clusters was obtained from A{B — V). 
Because of the approach of tying relative ages within lim- 
ited metallicity ranges to the absolute age of selected ref- 
erence clusters, we will, for the remainder of this paper, 
call our method simply the j4M-method. 

We have employed in our analysis the BV and VI 
isochrones described in SW98 and Weiss fc Salaris (1999 ), 
assuming a cosmological helium mass fraction Y=0.230, 
an average a-element overabundance [a/Fe]=0.4, and a 
helium enrichement factor AY/AZ=3. For the {V — /)- 
colours we use a combination of the transformations by 
[Bessell et al. (1998 ) for the main sequence and [von Braun 



et al. (1998) for the RGB, as described in Weiss & Salaris 
(1999[). For further det ails, see the original papers (SW98, 



ship given by CG97, but allows us to directly compare our 
results with R99, who selected the same [Fe/H] values for 
their sample. In case of clusters without estimates from 
Rutledge et al. (19971), but with estimates on the ZW84 



[Weiss fc Salaris 1999 ). We briefly recall here that the 
ZAHB luminosities from our models, which set the zero 
point of our absolute ages, are in agreement with the 
ZAHB brightness levels obtained by Carretta et al. (2000[ ) 
from the subdwarf fitting method applied to a sample of 
Galactic globular clusters. They are therefore about 0.13 
mag brighter than the relationship used by R99, which 
results from an average of the subdwarf fitting distances 
together with other d istance determinations, discussed in 
[Carretta et al. (2000[) . 

Our procedure allows a reliable determination of the 
age distribution with respect to the absolute age of the 
template clusters, largely unaffected by uncertainties in 
the convection treatment, by the precise value of the clus- 
ter metallicity, and by the inclusion of element diffusion 
(see the discussion in SW97) . It is also fairly independent 
of current uncertainties in colour transformations. This we 
verified by repeating the procedure for the clusters with 
VI photometry using a different set of RGB colour trans- 



scale, ^ ve again transformed their [Fc/H] onto the CG97 
scale using the conversion formula given by CG97. 

Within each interval a reference cluster showing a HB 
well populated at the instability strip region and/or at 
its red side was selected; this allows a reliable determina- 
tion of the ZAHB level to be obtained, and the absolute 
age to be determined accurately from AV^, by comparing 
it with AV taken from the isochrone for the metallicity 
under consideration. The reference clusters are, respec- 
tively, M 15 (NGC 7078), M 3 (NGC 5272), NGC 6171 
and 47 Tuc (NGC 104). In our new analysis M 15 and 



form ations together with our isochrones ( Alonso et al 



1999 ); we find only modest age changes within ±0.6 Gyr. 



A more detailed discussion on the issue of colour transfor- 
mations will be presented in the next section. 

Two further tests have been performed in order to as- 
sess the reliability of our age determination: In the first 
one we have simply moved clusters close to the boundary 
of two contiguous groups (e.g., NGC 6254, NGC 3201, 
NGC 6397) to the adjacent one and redetermined their 
age with respect to the absolute age of the new template 
cluster. The resulting absolute ages agree within 0.3 Gyr 
with the ones obtained from the original group location. 



4 



M. Salaris & A. Weiss: Ages of a large sample of Galactic globular clusters 





Fig. 1. Comparison of the absolute ages (in Gyr) obtained 
using the AV values provided by R99 for their whole sam- 
ple (vertical axis), with the ages obtained with our AM- 
method (horizontal axis) . Errors are as shown in Table |l|, 
or as given by R99. The solid line represents the locus of 
equal ages; dashed line corresponds to the average differ- 
ence between the two sets of ages (0.3 Gyr). 



In the second test, we have derived the age of all 35 clus- 
ters from R99 making use of their estimated AV^-values 
for the whole sample. Even in case of blue HB morpholo- 
gies, R99 have derived the level of the ZAHB extrapo- 
lating the observed ZAHB location towards the red, by 
fitting template cluster HBs uniformly populated to the 
observed branches. While we consider this approach be- 
ing less accurate then in cases with observed HBs pop- 
ulated at the instability strip or at its red side, it is of 
course a valid estimate of the ZAHB brightness. The re- 
sulting ages for the 35 clusters, displayed in Fig. |l|, indeed 
agree extremely well with the values obtained with the 
AM-technique, with an average difference of only 0.3 Gyr 
{AV ages being younger). This test confirms the mutual 
consistency of the AV relative ages and those obtained by 
the AM- method over the whole [Fe/H] range spanned by 
the cluster sample. 

A similar test, but in {B — V), was performed for two 
individual clusters, NGC 6712 and NGC 6934. BV pho- 
tometry is available for both of them, and both have ex- 
tremely well-defined ZAHB levels. We could therefore de- 
rive their absolute age also directly from AV. As in the 
case of the R99 cluster sample with VI photometry, the 
results are consistent, within 0.4 Gyr, with the relative 
age obtained from their A{B — V) with respect to the 
corresponding template cluster. 

The final results for the ages of the entire sample are 
displayed in Table ^ In case of non-template clusters, the 
final error in the absolute age has been derived by adding 



Fig. 2. Comparison of the absolute ages (in Gyr) ob- 
tained by SW98 (vertical axis) with the values found in 
the present work (horizontal axis) , in case of 24 common 
clusters for which we have used new photometric data. 
The solid line represents the locus of equal ages. 



in quadrature the error in the age of the corresponding 
template cluster (which is determined by the errors in AV 
and the internal errors in the adopted metallicity scale) to 
the error in the relative age estimate (which is determined 
by the errors in the differences in A{V — I) or A{B — V), 
and the internal errors in the adopted metallicity scale). 

There are 24 clusters in common with the SW98 sam- 
ple, for which we have now employed more recent and 
more homogeneous photometry (mainly in the VI bands); 
21 of these clusters belong to the R99 sample. In Fig. || 
we compare the absolute ages derived by SW98 (where 
essentially the same [Fe/H] scale was employed) and the 
ones we now obtained for these 24 clusters in common. 
There is a clear average shift towards higher values in our 
new sample, which is due to the new photometric data; 
the average difference is of the order of 1 Gyr. However, 
in case of two metal-rich clusters, namely NGC 6366 and 
NGC 1261 (new VI photometries from R99), we get ages 
considerably lower by ~2.5 Gyr. 

Fig. H displays the absolute ages as a function of [Fe/H] 
and galactocentric distance (Rgc) for our whole sample, as 
reported in Table |l|. Filled circles correspond to clusters in 
the R99 sample, while open squares denote the remaining 
ones. There are no significant differences between these 
two sub-samples (see also Sect. 4). In case of the clusters' 
Rgc we have adopted, as a reference, the values provided 
by Harris (1996). Strictly speaking, Rgc depends on the 
cluster distance modulus, but for differences of the order 



of 0.1-0.15 mag around the Harris (1996) values, Rgc is 
not changed appreciably, and in case of, e.g., the 35 clus- 
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Table 1. Data for the 55 clusters analyzed. The columns display, respectively, cluster name, age using the CG97 scale 
and associated error, [Fe/H] on the CG97 scale, age using the ZW84 scale and associated error, [Fe/H] on the ZW84 
scale, galactocentric distance, HB type, photometric bands used for the age determination, source of the data (see text 
for details). Reference clusters are emphasized by bold type characters. 



Cluster AgecG97(Gyr) [Fe/H]cG97 Agezw84(Gyr) [Fe/H]zw84 Rgc(Kpc) HBtypc phot source 



10.7 ± 1.0 
11.9 ± 1.1 
9.5 ± 
1.1 ± 



NGC 104 

NGC 288 
NGC 362 
NGC 1261 
NGC 1851 
NGC 1904 
NGC 2298 
NGC 2419 
NGC 2808 
NGC 3201 
NGC 4590 
NGC 5053 
NGC 5272 
NGC 5466 
NGC 5897 
NGC 5904 
NGC 6093 
NGC 6101 
NGC 6121 
NGC 6171 
NGC 6205 
NGC 6218 
NGC 6254 
NGC 6341 
NGC 6352 
NGC 6362 
NGC 6366 
NGC 6397 
NGC 6426 
NGC 6535 
NGC 6584 
NGC 6624 
NGC 6637 
NGC 6652 
NGC 6656 
NGC 6681 
NGC 6712 
NGC 6723 
NGC 6752 
NGC 6779 
NGC 6809 
NGC 6838 
NGC 6934 
NGC 7078 
NGC 7099 
NGC 7492 
Arp 2 
Eridanus 
IC 4499 
Pal 3 
Pal 4 
Pal 5 
Pal 12 
Rup 106 
Terzan 7 



10.7 
11.3 
8.7 
8.6 
9.2 
11.7 
12.6 
12.3 
9.3 
11.3 
11.2 



1.0 
1.1 
1.5 
1.1 
1.1 
± 1.3 
± 1.4 
± 1.0 
± 1.1 
± 1.1 
± 0.9 



10.8 ± 0.9 



11.3 
12.2 
12.3 
10.9 
12.4 
10.7 
11.7 
11.7 
11.9 
12.5 
11.8 
12.3 

9.9 
11.0 

9.6 



± 0.7 



0.9 
1.2 
1.1 
1.1 
1.4 
1.1 
± 0.8 
± 1.1 
1.3 
1.1 
0.9 
1.4 
1.3 
1.4 
1.1 



± 
± 
± 
± 
± 
± 

12.1 ± 
12.9 ± 1.0 
12.8 ± 1.2 
11.3 ± 1.4 
10.6 ± 1.4 
10.6 ± 1.4 

± 
± 
± 
± 

± 1.3 

12.2 ± 1.1 

12.3 ± 1.4 
12.3 ± 1.7 
10.2 ± 1.4 

± 
± 
± 
± 



11.4 
12.3 
11.5 
10.4 
11.6 



1.0 
1.2 
1.4 
1.4 



9.6 
11.7 
11.9 
12.0 
11.3 

8.9 ± 
12.1 ± 

9.7 

9.5 

9.8 

6.4 



1.5 
0.8 
1.4 
1.4 
± 1.4 
1.6 
1.4 
1.3 
1.6 
1.4 
0.9 



± 
± 
± 
± 

10.2 ± 1.4 
7.4 ± 1.4 



-0.78 
-1.14 
-1.09 
-1.08 
-1.03 
-1.37 
-1.71 
-2.14 
-1.11 
-1.24 
-2.00 
-1.98 
-1.33 
-2.13 
-1.73 
-1.12 
-1.47 
-1.76 
-1.05 
-0.95 
-1.33 
-1.14 
-1.25 
-2.10 
-0.70 
-0.99 
-0.73 
-1.76 
-2.11 
-1.51 
-1.30 
-0.70 
-0.78 
-0.81 
-1.41 
-1.35 
-0.94 
-0.96 
-1.24 
-1.61 
-1.54 
-0.73 
-1.30 
-2.02 
-1.92 
-1.41 
-1.45 
-1.20 
-1.26 
-1.39 
-1.07 
-1.24 
-0.83 
-1.49 
-0.56 



9. 



1.5 
1.1 



9.1 ± 1.1 
12.6 ± 1.3 



12.9 
12.8 
10.2 
12.1 



± 1.4 

± 1.0 

± 1.1 

± 1.1 



11.2 ± 0.9 



10.1 
12. 



0.9 
0.7 



± 
± 

± 0.9 
± 1.2 
± 1.1 
± 1.1 

± 
± 

11.7 ± 0.1 

± 
± 
± 
± 
■ ± 
± 



12.5 
12.4 
11.6 
12.9 
11.0 
11.9 



13.0 
12.7 
12.2 
12.8 
9.7 
11.1 



1.4 
1.1 



1.3 
1.3 
1.1 
0.9 
1.4 



9.4 ± 
12.5 ± 



13.0 
13.1 



12.1 ± 
10.6 ± 



1.3 
1.4 
1.1 
1.0 
1.2 
1.6 



1.4 



10.6 ± 

11.4 ± 

12.5 ± 
11.9 ± 

10.5 ± 

11.6 ± 

12.7 ± 

12.8 ± 
12.4 ± 
10.1 ± 
10.0 ± 
11.8 ± 
12.3 
12.1 



1.4 
1.0 
1.2 
1.4 
1.4 
1.3 
1.1 
1.4 
1.7 
1.4 
1.6 



± 1.4 
± 1.4 
11.5 ± 1.4 
8.4 ± 1.6 
11.2 
9.2 
9.2 ± 1.6 
10.0 
6.4 
10.4 
7.5 



± 1.2 
± 1.3 



± 1.4 
± 0.9 
± 1.4 
± 1.4 



-0.71 
-1.40 
-1.33 
-1.32 
-1.23 
-1.67 
-1.85 
-2.10 
-1.36 
-1.53 
-2.11 
-2.10 
-1.66 
-2.22 
-1.93 
-1.38 
-1.75 
-1.81 
-1.27 
-1.09 
-1.63 
-1.40 
-1.55 
-2.24 
-0.50 
-1.08 
-0.58 
-1.94 
-2.20 
-1.78 
-1.54 
-0.50 
-0.72 
-0.89 
-1.75 
-1.51 
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-0.58 
-1.54 
-2.13 
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-1.47 
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-1.90 
-1.00 
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0.52 
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0.86 
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0.93 
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0.97 
0.97 
0.98 
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-0.58 
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0.98 
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0.91 
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-0.08 
1.00 
0.98 
0.87 
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0.25 
0.67 
0.89 
0.81 
0.86 
-1.00 
0.11 
-0.50 
-1.00 
-0.40 
-1.00 
-0.82 
-1.00 
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Testa et al. (200li TOl), [Harris et al. (1997|, H97), ISarajedini fc Da Costa (1991 



H99), ^arajcdini fc Forrester (199E|, SF95) , |Heasley et al. (200C|, HOO), [Chaboyer et al 



SDC91), [Hatzidimitriou et al. ( 1999 
2000|, COO), iPaltrinieri et al (2001 



POl), iPiotto et al. (19991, P99), [Richer et al (198^ , RFV88), [Cote et al. (1991| , CRF91), [B uonanno et al. (1998| , B98), ^tetson 
et al. (199£| , S99), ^mith et al. (198e| , SM86). 
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count the error coming from the use of A(y — /) as the 
relative age indicator, without adding the contribution of 
the error in the absolute age of the template cluster (as is 
done in Table |^), because for typical absolute age errors 
of ±1 Gyr, the differential properties of A(y — /) as a 
function of metallicity and age are negligibly affected. 

For the pair NGC 288-NGC 362 we obtain an age dif- 
ference of 2.6±1.5 Gyr (NGC 288 being older), in good 
agreement with the results by Bellazzini et al. (2001), but 
in contrast to SW97, where the two clusters appeared to 
be coeval. The reason for this change lies in the new pho- 
tometric data. According to the simulations by patelan 



et al. (2001), this age difference, coupled with the abso- 



lute ages we obtain and the use of the CG97 [Fe/H] scale, 
can explain the overall different HB morphology; however, 
canonical HB models appear unable to reproduce the de- 
tailed morp hology of the red end of NGC 288 HB ( |C at elan 
et al., 200l| ). 



log(Rgc(Kpc)) 

Fig. 3. Distribution of the absolute ages (in Gyr) as func- 
tion of [Fe/H] (upper panel) and Rgc (lower panel). Filled 
circles correspond to clusters in the R99 sample, while 
open squares denote the remaining ones. 




0.5 



1.0 1.5 
log(Rgc(Kpc)) 



2.0 



Fig. 4. As in Fig. |, but assuming the ZW84 [Fe/H] scale 
for the clusters instead of the CG97 one in the previous 
figure. 



ters with ZAHB levels from R99, we determined distance 



moduli within this range from the Harris (1996) values. 

It is interesting to determine the age difference for 
two classic second-parameter pairs, namely NGC 288- 
NGC 362 and M 3-M 13 (NGC 6205), both contained in 
our sample. When deriving errors in age differences within 
the same metallicity group, we have just to take into ac- 



In case of M 3-M 13 we obtain a difference of 0.6±1.0 
Gyr (M13 older), much less signific ant than the 1.7±0.7 
Gyr as obtained by Rey et al. (200l| ) from their BV pho- 
tometry. M 13 was not contained in our previous work. 
The same negligible age difference we obtain from the 
Ay values derived by R99; it seems therefore that the 
discrepancy between our and Rey et al. (2001) result is 
due to real differences in the photometric data, and not 
to the use of different passbands for the TO-RGB colour 
differences and inconsistency in the colours of the theo- 
retical isochrones. Problems with the calibration of the 
photometry may possibly lead to this kind of discrepancy. 



Rey et al. (2001) noticed, for example, that their derived 



fiducial line for M 13 agrees well with those obtained by 
Richer fc Fahlman (1986| ) and [Yim et al. (2000| ), but dif- 
fers in the mai n sequence and subgiant branch region from 
the fiducial by Paltrinicri et al. (1998 ). 

It is also interesting to notice that Arp 2 and Rup 106 
do not appear much younger than the bulk of the clusters 
at their metallicity; the reason why Buonanno et al. (1998 ) 
found higher age differences is mainly the fact that their 
adopted clusters ' absolute ages are higher. As discussed 
in detail in, e.g., Pulone et al. (199^ ), lower absolute ages 
imply smaller age differences for a given observed distribu- 
tion of Ay, A(i? — V) or A(y — /) values. As discussed in 
SW98 there are preliminary indications that Rup 106 and 
also Pal 12 may not show a-element enhancement. In this 
case, their ages displayed in Table should be increased 
by about 1 Gyr. 

To highlight the effect of the present uncertainties in 
the [Fe/H] scale, we have also derived ages, as a test, by 
using the [Fe/H] values given by Rutledge et al. (1997) on 
the ZW84 scale (internal accuracy again of the order of 
0.10 dex), complemented, if needed, by data in the orig- 
inal ZW84 paper or coming from other spectral indices 
calibrated on the ZW84 scale. In this case the [Fe/H] 
range spanned by our sample is larger than the CG97 
scale. We have therefore divided the sample into 5 groups, 
having as template clusters M 15 (-2.3 <[Fe/H]< -2.0), 
NGC 6656 (-1.99 <[Fc/H]< -1.7), M 3 (-1.69 <[Fe/H]< 
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-1.4), NGC 6171 (-1.39 <[Fe/H]< -0.9) and 47 Tuc 
(-0.89 <[Fe/H]< -0.5). Fig. |shows the age distribution 
resulting from the AM- method as a function of [Fe/H] and 
Rgc, when using the ZW84 [Fe/H] scale. Using this metal- 
licity scale, the age differences between the two second 
parameter pairs discussed remain basically unchanged. In 
fact, in case of NGC 288-NGC 362 we obtain now a dif- 
ference of 2.4±1.5 Gyr, and for M 3-M 13 we get 0.9±1.0 
Gyr. A detailed analysis of the age distributions displayed 
in Fig. ^ and Fig. ^ will be presented in the last section. 

3. Colour-based relative ages 

We consider it important that isochrones - together with 
colour transformations - yield ages largely independent 
of the age indicator. If consistency is achieved, one could 
then in principle apply different age indicators to different 
clusters, according to the quality of the available data, and 
still have consistent ages for all the sample. 

In the last section we showed that the AM-method 
yields results which are consistent with absolute and rel- 
ative ages from the AV - or vertical method. Another 
age indicator - the horizontal one - is the colour differ- 
ence between TO and RGB, which we used already in the 
AM-method to obtain relative ages within each metallic- 
ity bin. In the following, we will call this quantity simply 
SC, independent of the actual colour it refers to. From an 
observational point of view this age-indicator is the easi- 
est one to determine, since it uses a CMD portion which 
is common to all clusters, is well-populated, and thus - 
for sufficiently deep photometry - observable with high 
accuracy. Reddening is not important, as long as it is low 
and not differential. However, the difficulties start when 
using it for absolute ages or over a wide range of metallic- 
ities; in this case, observed individual 6C values have to 
be compared to theoretically predicted ones, and colour 
transformations ~ with their known uncertainties - must 
be able to predict with extreme accuracy the absolute val- 
ues of colours over a large [Fe/H] and effective temperature 
range. 

Because of the eas y observational accessibility, 
Buonanno et al. (199§| ) developed a semi-empirical 



relative-age determination method using 6C in the BV 
plane. It sets out by selecting a group of coeval clusters 
of all metallicities, which defines the empirical function 
^C([Fc/H]) at the age to of this group. For this step, 
AV and theoretical isochrones are needed. It turns out 
that membership to the coeval group is rather indepen- 
dent of the isochrones used, but not the actual value for 
tg. The result for our own isochrones is shown in Fig. ^. 
The empirical (5C( [Fe/H]) -relation at can be compared 
to the theoretical one. Agreement is not necessarily given 
(Buonanno et al., 1998) and depends on isochrone set and 




Fig. 5. Observed AV values for the 35 clusters from R99 
along with the theoretical values from our isochrones, plot- 
ted for 14-9 Gyrs (bottom to top in steps of -1 Gyr). 



obtains SC as a function of age and metallicity, such that 
the observed values for SC along with [Fe/H] yield the age 
straightforwardly. This makes the method very attractive 
for large cluster samples. Because of the problems men- 
tioned above the method works best for relative ages, but 
one has to keep in mind that age differences also depend 
on the absolute age (SW97). In [Pulone et al. (1998| ) we 
showed that our own isochrones and colour transforma- 
tions to {B — V) result in a (5C([Fe/H]) relation consistent 
with observations and that relative ages from this method 
agree well with those we obtain from the AM-method, 
which mi ght be understood as a "piecewise" application 
of that of 



Buonanno et al. 



This same idea has then been used by R99 for their 
clusters observed in V and / (Rosenberg et al., 2000a 
They employed isoc hrones from ^traniero et al. (1997| ) 
and a pre-release of V^andenBerg et al. (2000| ) models for 
the theoretical determination of 6C{t) at various [Fe/H], 
and the TO brightness at various ages and [Fe/H]; as for 
the HB brightness, they considered the empirical ZAHB 



brightness relationship by Carretta et al. (2000), lacking 
access to ZAHB models consistent with the isochrones 
adopted. 



colour transformation. 



As in Buonanno et al. (1998) the comparison of ages 
from AV and A(V^ — /) revealed inconsistencies: the ab- 
solute age to of the coeval reference sample depends on 
age indicator and isochrones: it is (ages in Gyr) ^^14 for 
VandenBerg et al. (2000) and ^15 for Straniero et al 



In ihe next step clusters of ages differing from to are 



used to derive the relation dC{t) at various [Fe/H] points. 
Age differences are again determined from the vertical 
method and isochrones. Combined with tlie first step, one 



(1997), according to AV^, but, respectively, 16 and 13 for 
A(V^ — /), i.e. there are changes of 2 Gyr with respect to 
the AV ages, but they are in opposite directions for the 
two isochrone sets. 
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Fig. 6. The observed horizontal relative age indicator 
A{V — I) for the 35 clusters from R99 along with the the- 
oretical values from our isochrones, plotted for 14~9 Gyrs 
(bottoin to top in steps of -1 Gyr). Isochrones were trans- 
formed to {V ^I) using th e colour transformation selected 



in [Weiss fc Salaris (lyijgU 



We already demonstrated (Fig. |l|) that the relative 
AM-ages, which are basically obtained from the horizon- 
tal method, agree well with the AV^ values for the R99 
sample. In particular, the majority of clusters is coeval at 
« 12 Gyr with a small metal-rich group at ages lower by 
1-3 Gyr (see Fig. |^). According to the (global) horizontal 
age indicator A{V — /), however, there is a clear age in- 
crease of almost 3 Gyr within the coeval group (Fig. 
This effect is most pronounced at the lowest metallicities. 
The A{V — /) therefore, when used as an age indicator 
across all the metallicity range, along with our isochrones 
and colour transformations, produces inconsistent results 
with respect to AV. Note, however, that the horizontal 
average age of the "coeval" group is around 11 Gyr and - 
within the errors - consistent with that of the vertical one 
(Fig. H), and that the whole discrepancy in the A{V — I) 
relative ages is due to a variation in A{V — I) of only 
« 0.05 mag. This is a strong indication that theoretical 
colours have to be very accurate for the whole metallicity 
range of globular clusters in order that this method can 
return reliable results. 

To emphasize this point we show in Fig. 7 the predicted 
RGB colours for 12 Gyr-isochrones of different metallici- 



ties as obtained by using the transformations of Weiss & 
Safaris ] (1999| ) or [Alonso et al. (1999|). Evidently, in th e 
former case, which is based on von Braun et al. (1998| ) 



the dependence of colour on metallicity is non-monotonic 
and RGBs are piling up at a certain bluest level, while 



for Alonso et al. (1999) the behaviour is much more regu- 
lar. Taking this as the "true" behaviour, von Braun et al. 
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Comparison of RGB {V — I) colours for isochrones 
G[yr and [Fe/H]=-2.3, -2.0, -1.6, -1.3, -0.7 (from 
ht). Solid lines correspond to the colour trans- 
m selected in |Weiss & Salar is (1999 ), dashed ones 



by Alonso et al. (199£) 



(1998) RGBs would be too red for the lowest metallicities 
by about 0.04-0.06 mag, i.e. ages would be too high by 
the corresponding age difference of 2-3 Gyr, just as vis- 
ible in Fig. p. Note that at higher metallicities the two 
transformations agree well. Repeating therefore the hor- 
izontal absol ute age determinatio n, but using the trans- 
formation by Alonso et al. (1999| ), we obtain the results 



displayed in Fig. y. While there is still a small age gradient 
across the main group of clusters, the total age difference 
is reduced to within 1.5 Gyr, and is consistent within the 
errors with the results of the AM-mcthod - which, as al- 
ready discussed in the previous section, are unaffected by 
the colour transformation choice - and from the AV. This 
shows that our isochrones together with the colour trans- 
formations by Alonso et al. (199^ ) for the RGB are able 
to produce stable relative ages for the whole metallicity 
range of GCs independent of the use of the AM-method, 
the AV^, or 6C age indicator. 

4. Discussion 

Having fully assessed the consistency of the age distri- 
bution within our cluster sample, the results displayed in 
Fig. ^ can be used to investigate the formation mechanism 
of the Galaxy. In the following we update the discussion 
in SW97 and SW98, taking into account this much larger 
cluster sample. 

There are about 150 globular clusters known in the 



Galaxy (Harris 1996), and probably the total population 
is of at most about 200 objects (Harris 1976); therefore, 
our sample contains between 25% and 35% of the Galactic 
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Fig. 8. As Fig. b ut with the colour transformations by 
Alonso et al. (1999[ ). 



globular cluster population, and should provide statisti- 
cally more significant data about the Galaxy formation 
mechanism, with respect to previous works. 

Two fundamental pieces of information need to be ex- 
tracted from our age-analysis, namely the age distribution 
of the clusters with respect to their initial chemical com- 
position (whose diagnostic is [Fe/H]), and with respect 
to their position within the Galaxy (whose diagnostic is 
Rgc). To detect an intrinsic age spread - that is, not due to 
the error bar on the age determinations - among a given 
sample of objects which shows an age dispersion (Jobs, we 
have followed the same procedure as in SW97 (see also 



Chaboyer et al. 1996). More in detail, we calculate an "ex- 



pected" distribution for the assumption of no intrinsic age 
range (crrangc=0) by randomly generating, for each object, 
10000 ages using a Gaussian distribution. The mean of the 
distribution is given by the mean age of the clusters in the 
given sample, and the Ccxp by the typical error in the in- 
dividual age determinations (in this case we consider the 
total error given in Table ^ since we compare ages across 
multiple metallicity ranges). This is repeated for all clus- 
ters in the sample, and the F-test is then applied to the 
observed age distribution and this synthetic one, to test 
if they both show the same variance. We state that an 
age range exists if the probability that the two distribu- 
tions have the same variance is smaller than 5%. If this 
is true, the intrinsic age dispersion is then estimated by 



2 

obs 



^0.5 



4.1. Age-metallicity relationship 

When considering the complete sample of 55 objects, we 
obtain that all clusters with [Fe/H]< —1.2 result to be 
coeval (within the stated error bars), with an average age 



(t)=11.7±0.7 Gyr, where the dispersion of 0.7 Gyr is en- 
tirely due to the errors in the age determination. In case 
of [Fe/II]> —1.2 a statistically significant age spread is 
present. The average age is (t)=10.1 Gyr, with an intrinsic 
spread crrangc=l-9 Gyr. In this same [Fe/H] range we have 
detected, overimposed to the age spread, a marginally sig- 
nificant age gradient ^t/^[Fe/H]=— 2.4±1.8 Gyr/dex. 

Our complete sample contains 3 clusters generally 
considered to be part of the thick disk of the Galaxy 
(NGC 104, NGC 6352, NGC 6838), and 2 clusters pos- 
sibly belonging to the bulge population (NGC 6637, 
NGC 6324). All of these five clusters are metal rich and lo- 
cated in the inner halo. If we repeat our previous analysis 
excluding these 5 objects, and therefore considering only 
bona fide halo clusters, we obtain an average age (t}=10.0 
Gyr for the metal rich group, with crrangc=2.5 Gyr, and a 
more significant Jt/J[Fe/H]=— 4.3±2.4 Gyr/dex. 

4.2. Age-galactocentric distance relationship 

The clusters in the so-called inner halo (Rgc < 8 kpc) ap- 
pear to be coeval, with (t) = 11.4±0.8 Gyr, while the outer 
halo clusters display an intrinsic age range (Tiango=2.5 Gyr. 
The average age of the outer halo sample is (t) = 10.7 Gyr, 
lower than the inner halo one, simply due to the age spread 
towards lower ages, but there is no significant age gra- 
dient in the outer halo. This would imply that globular 
clusters have started forming at the same time through- 
out the Galaxy - this holds of course in case the actual 
Rgc is a good approximation of the value at the beginning 



of Galactic formation; see, e.g., the discussion by Searle 
|fc Zinn (1978 ) on this subject - with younger and more 
metal rich objects being formed or possibly accreted only 
at large Rgc values. This scenario is very similar to the one 
proposed by Searle fc Zinn (197S ), where the inner halo 



collapsed on short timescales, while the outer halo under- 
went a much longer formation period, accreting material 
over several Gyr. 

When considering only bona fide halo clusters the in- 
ner halo appears still to be coeval, its average age being 
basically unchanged ((t)=11.7±0.6 Gyr). 

4.3. Influence of the sample size and metallicity scale 

The previously described results about the age distribu- 
tion with respect to [Fe/H] and Rgc are unchanged if we 
restrict our analysis to the more homogeneous sample of 
35 clusters from R99. Moreover, the age trends with re- 
spect to [Fe/H] and Rgc we derive are similar to the re- 
sults obtained by R99 with their different technique, and 
statistically more significant because of the larger clus- 
ter sample. When comparing our new results with SW98 
we find that the general picture of age versus [Fe/H] is 
not greatly modified. In our previous works we obtained 
an approximate constant age for the oldest clusters at all 
metallicities, plus the onset of an age spread above a given 
[Fe/H]. However, while in SW97 and SW98 we did not find 
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any clear correlation between age and Rgc, we find now a 
difference between the age distribution of the inner and 
outer halo clusters. This is due to the much larger clus- 
ter sample and also to the upward revision of the age of 
NGC6652, which in SW97 and SW98 resulted to be much 
younger than the other inner halo clusters. 

The details of the scenario emerging from our age 
determination are however dependent on the selected 
[Fe/H] scale. When we consider the ages obtained with 
the ZW84 scale displayed in Fig. ^ we find that only 
clusters with [Fe/H]< —1.6 are coeval, with an average 
age (t)=12.2±0.6 Gyr, slightly larger than in the case of 
the CG97 metallicities. For larger [Fe/H] values an intrin- 
sic spread o'rangc=2.1 Gyr is present, with an average age 
(t)=10.5 Gyr, and (5t/(5[Fc/H]=-1.3±0.7 Gyr/dex. Also 
the distribution of ages as a function of Rgc is more homo- 
geneous, with an intrinsic age spread appearing only when 
Rgc is larger than ~13 kpc. If we restrict the analysis to 
only bona fide halo clusters, we find for [Fe/H]> —1.6 a 
slope 5t/5[Fe/H]=— 2.3±1.1 Gyr/dex, the significance of 
which now exceeds the 2a level. 

These differences come obviously from the different 
metallicity distribution of the clusters. In particular, in 
the range —1.8 <[Fe/H]< —1.0, large differences (up to 
'^O.S dex) exist between the two scales. 

Taking into account the effect of the uncertainty in 
the [Fe/H] scale and the cluster sample size, helps in ex- 
plaining discrepant conclusions reached by different au- 
thors about the existence and the onset of the age spread 
in the cluster population. As a test, we have considered 
the clusters used by VandenBerg (20do| , VOO) in his re- 
cent analysis; by applying our statistical analysis to his 
age determinations, wc find the onset of an age variation 
among clusters starting at [Fe/H]^ —1.6, lower than the 
value [Fc/H]=-1.2 found by R99, who used the CG97 
scale, and also lower than our results when using the CG97 
scale. A significant age-metallicity relationship is found for 
[Fe/H]> -1.6. 

VOO derived ages for 26 objects - 25 of which are also 
included in our sample - by employing cluster [Fe/H] val- 
ues close to the ZW84 scale, and absolute age determi- 
nations based essentially on At^ (more precisely, the dis- 
tance modulus is determined from fitting the appropriate 
theoretical ZAHB to the observational counterpart, and 
then the age is determined by the isochronc which best 
fits all the TO region of the observed Colour-Magnitude- 
Diagram) . 

We have rederived the ages for the 25 clusters in com- 
mon assuming VOO [Fe/H] values, our isochrones, our 
method and our adopted photometries (which are dif- 
ferent from the data employed by VOO for the majority 
of the clusters), and the results are displayed in Fig. ||, 
together with VOO ages shifted by —1 Gyr. When this 
shift of the absolute age is applied ~ which rcficcts mainly 
the ZAHB brightness difference between the different sets 
of models employed - the age distribution appears quite 
similar, all pairs of values being in agreement within 
the respective error bar. This consistency is also con- 




Fig. 9. Comparison between the ages derived with the 
AM-method, our isochrones and our adopted photome- 
tries (filled circles), and the ages from VOO (open squares) 
for 25 clusters in common; an offset by —1 Gyr has been 
applied to the VOO data. For the sake of clarity VOO error 
bars (of the order of 0.8 Gyr) have been omitted. 



firmed by the application of a K-S test to the two sam- 
ples. If we restrict the analysis to bona fide halo clus- 
ters (thus not considering 47 Tuc in the VOO sample), for 
[Fe/H]> —1.6 we obtain an age-metallicity relationship 
with slope (5</(5[Fe/H]=— 4.7±0.8 Gyr/dex, in good agree- 
ment with the value (5t/(5[Fe/H]=— 5.7±1.2 Gyr/dex one 
obtains from VOO ages. This slope is about a factor of 
two higher than the one we obtained with our full sample 
of bona fide halo clusters coupled with the ZW84 [Fe/H] 
scale. 

4.4. Second parameter and age 

With our large cluster sample it is also possible to investi- 
gate the global correlation between [Fe/H], HB type and 
age for the Galactic globular cluster system as a whole. 
We have already discussed in Sec. 2 two specific pairs of 
second-parameter clusters; in case of NGC 288-NGC 362 
we found that age could be the main reason for the dif- 
ferent HB types, while we do not find a significant age 
difference for the pair M 3-M 13. 

In order to detect possible correlations between age 
and HB type (defined as HBtypo=(B-V)/(B-|-V-|-R), where 
B,V,R are the number of HB stars located, respectively, 
at the blue side of the instability strip, in the instability 
strip and at its red side) superimposed on the relation- 
ship between HB type and [Fe/H], we have performed the 
following test. We have first considered the samples of co- 
eval clusters with Rgc < 8 for the CG97 [Fe/H] scale, and 
Rgc < 13 for the ZW84 one; since in both cases they span 
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Fig. 10. Relationship between age differences A(t) and 
normalized HB type for clusters in the [Fe/H] ranges (on 
the CG97 scale) indicated in the three panels. See text for 
details. 



4 
2 

-2 

-4 - 
4 

2 



-2 

-4 - 
4 

2 



-2 

-4 - 



o 



-1.6 < [Fe/HU, S-O.g 

gO o o 
o 



o 
o 



[Fe/H]„ 



-0.9 



-0.5 



0.5 



Fig. 11. As in Fig. ^ but considering the ZW84 [Fe/H] 
scale. See text for details 



the entire metallicity range covered by the whole sample, 



we use thcni to determine the variation of HBtypc with 
respect to [Fe/H], at constant age. The entire sample has 
then been divided into different metallicity bins, specified 
in Fig. |l^ and Fig. |ll|. For the most metal poor bin the 
average age of the sample has been subtracted from the 
individual cluster ages; in case of the two most metal rich 
bins, the average age of their combined sample has been 
subtracted from the individual values. These normalized 



individual ages are denoted as A(t). To remove the in- 
fluence of the first parameter ([Fe/H]), in addition, a ref- 
erence [Fe/H] has been selected for each range, and the 



individual values of HBtypo (taken from Harris 1996 and 



displayed in Table |l]) have been corrected using the deriva- 
tive AHBtypc/A [Fe/H] (determined from the previously 
mentioned coeval samples) and considering the difference 
between the actual cluster [Fe/H] value and the refer- 
ence one. By plotting A{t) as function of HB^yp^,, that is, 
the observed HB type normalized to the reference [Fe/H] 
value, one should in principle be able to isolate an hypo- 
thetical correlation between age differences and HB colour 
differences, independent of the effect of [Fe/H] variations 
among clusters. The reference [Fe/H] values are, respec- 
tively, [Fe/H]=-1.75 and -1.1 for the CG97 scale, and 
[Fe/H]=-2.0 and -1.3 for the ZW84 one. We did not cor- 
rect the HB type for the most metal rich cluster groups 
in Fig. |l^ and Fig. [l|, since the value of HBtypc saturates 
and does not provide useful information about a possible 
correlation with age. 

The upper panels of both Fig. |l^ and Fig. |l| display 
[Fe/H] ranges where our F-test analysis did not detect 
a significant age spread compared to the errors in the 
age estimates. In these ranges there is some spread in HB 
type but no significant correlation with age is found. In 
the middle panels the F-test analysis detected significant 
age spreads and also a statistically significant correlation 
with HB type is found. The bottom panels display clusters 
showing an age spread but their in this [Fe/H] range the 
HBtypc parameter saturates and no relevant information 
can be deduced. 



5. Summary 

We have determined absolute and relative ages for a large 
sample of globular clusters, using new and homogeneous 
photometric data, mostly in V and /. We show that dif- 
ferent age determination methods result in very similar 
ages, if colour transformations are carefully chosen for 
those methods using colour differences. Our preferred AM- 
method appears however insensitive to the colour trans- 
formations choice. In a few cases, new photometric data 
have led to a substantial change in age compared to our 
previous works; this source of uncertainty should not be 
ignored. 

From the analysis of our large sample of 55 Galactic 
globular clusters, we have found that: 

i) Irrespective of the adopted [Fe/H] scale, the metal 
poorer clusters are coeval within '^1 Gyr, while an age 
spread appears at higher metallicites. The onset of the 
age spread is at [Fe/H]'--^ —1.2 when using [Fe/H] values 
on the CG97 scale, and [Fe/H]~ —1.6 if one adopts the 
ZW84 scale. 

ii) An age-metallicity relationship exists among the 
non-coeval clusters. In case of the CG97 scale we 
find (5t/(5[Fe/H]=— 2.4±1.8 Gyr/dex when considering the 
complete sample; if we restrict the analysis only to bona 
fide halo clusters we obtain a more significant slope 
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5t/(5[Fe/H]=-4.3±2.4 Gyr/dex. In case of the ZW84 scale 
we obtain Jt/J[Fe/H]=— 1.3±0.7 Gyr/dex for the complete 
sample, and St /S[Fe/B.]=— 2. 3±1.1 Gyr/dex for the pure 
halo sample. 

iii) Irrespective of the adopted [Fe/H] scale, clusters 

closer to the galactic centre arc coeval, while an age spread 
appears at larger galactocentric distances. In case of the 
CG97 [Fe/H] scale the age spread starts at Rgc larger than 
^ 8 kpc, while in case of the ZW84 scale it starts at Rgc 
above 13 kpc. In both cases there is no significant age- Rgc 
relationship among clusters showing an age spread. 

iv) Our results yield indications that age differences 
are one main factor able to explain the global behaviour 
of the second parameter effect. 
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